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Abstract 

We present theoretical studies of high-order harmonic generation (HHG) produced by non- 
homogeneous fields as resulting from the illumination of plasmonic nanostructures with a short 
laser pulse. We show that both the inhomogeneity of the local fields and the confinement of the 
electron movement play an important role in the HHG process and lead to the generation of even 
harmonics and a significantly increased cutoff, more pronounced for the longer wavelengths cases 
studied. In order to understand and characterize the new HHG features we employ two differ- 

C/3 ■ 

t>-j" ent approaches: the numerical solution of the time dependent Schrodinger equation (TDSE) and 

in " 

Q-f the semiclassical approach known as Strong Field Approximation (SFA). Both approaches predict 

comparable results and show the new features, but using the semiclassical arguments behind the 
SFA and time-frequency analysis tools, we are able to fully understand the reasons of the cutoff 
extension. 
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I. INTRODUCTION 



Coherent light sources in the ultraviolet (UV) to extreme ultraviolet (XUV) spectral 
range are in high demand nowadays for basic research, material science, biology and possibly 
lithography l|. Their caveat is a demanding infrastructure for XUV generation and target 
delivery as well as its low efficiency and low duty cycle. The recent demonstration based 
on surface plasmon resonances as light enhancers could provide a possible solution to this 
problem 

Field enhanced high-order harmonic generation (HHG) using plasmonics, generated start- 
ing from engineered nanostructures, requires no extra cavities or laser pumping to amplify 
the pulse power. By exploiting surface plasmon resonances, local electric fields can be en- 
hanced by more than 20dB [sl,!^]. This amplification is strong enough to exceed the threshold 
laser intensity for HHG generation in noble gases and the pulse repetition rate remains un- 
altered without any extra pumping or cavity attachment. Furthermore, the high harmonics 
radiation generated from each nanostructure acts as a point-like source, enabling collimation 
or focusing of this coherent radiation by means of (constructive) interference. This opens a 
wide range of possibilities to spatially arrange nanostructures to enhance or shape spectral 
and spatial properties in numerous ways. 

The basic principle of high-order harmonic generation (HHG) based on plasmonics can be 
summarized as follows (the full explanation can be found in j2|): a femtosecond low intensity 
pulse is coupled to the plasmon mode inducing a collective oscillation of free charges within 
the metal. The free charges redistribute the electric field around each of the nanostructure, 
thereby forming a spot of highly enhanced electric field. The enhanced field exceeds the 
threshold of HHG, thus by injection of noble gases onto the spot of the enhanced field, high 
harmonics are generated. 

I„ the seminal expen,„e„t of K,. et al. g the output bea™ emitted fo. a modest power 
femtosecond oscillator (100-kW peak power, 1.3 nJ pulse energy, 10 fs pulse duration and 
800 nm of wavelength) was directly focused onto a 10 /imxlO fim bow-tie nanoantenna array 
with a pulse intensity of 10^^ W-cm~^, which is about two orders of magnitude less than 
the threshold intensity to generate HHG in noble gases. Their experimental result showed 
that the field intensity enhancement factor exceeded 20 dB, which is sufficient to produce 
XUV wavelengths from the 7th (114 nm) to the 21st (38 nm) harmonics with the injection 
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of xenon gas. Additionally, each bow-tie acts as a point-like source, thus 3D arrangement 
of bow-ties should enable us to perform control of generated harmonics in various ways by 
exploiting interference. 

Numerical and semiclassical approaches to study laser-matter processes in atoms and 
molecules, in particular High-order Harmonic Generation (HHG), are largely based on the 
assumption that the laser electric field is homogeneous in the region where the electron 
dynamics takes place Due to the strong confinement of plasmonic hot spots, the 

laser electric field is clearly not homogeneous anymore, in the region where the electron 
dynamics takes place, and consequently important changes in the laser matter processes 
would occur. From a theoretical viewpoint, the HHG process can be tackled using differ- 
ent approaches (for a summary see e.g. [3, [sl and references therein). In this article we 
concentrate our effort in extending two of the most and widely used approaches: the nu- 
merical solution of Time Dependent Schrodinger Equation (TDSE) in reduced dimensions 
and the Strong Field Approximation 9[. So far theoretical work in strong field physics has 
been focused on the processes driven by homogeneous coherent electromagnetic radiation 
on atoms and molecules. Only recently studies about how HHG spectra are modified due 
to non-homogeneous fields, as those present in the vicinity of a nanostructure irradiated by 



a short laser pulse, have been published jlOl. 



The paper is organized as follows. In the next two sections (Sec. II. A and Sec. II. B) 
we will present two theoretical approaches to model high-order harmonic generation (HHG) 
spectra produced by non-homogeneous fields, namely the numerical solution of the Time 
Dependent Schrodinger Equation in reduced dimensions (ID-TDSE) and the Strong Field 
Approximation (SFA), respectively. We have developed both approaches in such a way to 
allow the treatment of very general non- homogeneous fields showing the flexibility of both 
the ID-TDSE and SFA models. Particular studies of HHG spectra for the most simplest 
case is presented in Sec. HI. We will discuss how the non-homogeneous field produce no- 
ticeable modifications in the HHG spectra, namely a change in the harmonic periodicity 
(odd and even harmonics now appear) and a distinct extension in the position of the cutoff 
(more pronounced for the longer wavelengths cases studied). Both the ID-TDSE and SFA 
approaches provide comparable predictions, but analyzing the HHG process using semiclas- 
sical arguments we can present strong evidence about the motives of the cutoff extension. 
The paper ends with a short summary and an outlook. 
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II. THEORY 



A. Numerical solution of the Time Dependent Schrodinger Equation in reduced 
dimensions (ID-TDSE) 

Considering the dynamics of an atomic electron in a strong laser field is mainly along 
the direction of the field (in a linearly polarized laser pulse), it is reasonable to model the 
high-order harmonic generation (HHG) in a ID spatial dimension by solving the following 
Schrodinger equation (ID-TDSE): 

i^^^ = nm{x^t) (1) 
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To model an atom in ID, it is common to use the quasi-Coulomb potential 

Vatom{x) = ^== (2) 



which was first introduced in [12| and has been widely used in the ID analysis of an atom. 
The potential due to the laser electric field linearly polarized along the x-axis will be modified 
in order to treat nonhomogeneous fields. Consequently we write 

Vlaser{x,t) = E{x,t)x (3) 

with 

E{x,t) = Eof{t){l + eg{x)) smut (4) 

the laser electric field. In (jl]) Eq is the peak amplitude and u the frequency of the coher- 
ent electromagnetic radiation. Furthermore, f{t) defines the pulse envelope and e a small 
parameter that characterizes the inhomogeneity region. g{x) represents the functional form 
of the nonhomogeneous field. In this work we concentrate our analysis in the simplest form 
for g{x), i.e. g{x) = x, but we should to emphasize that the numerical algorithm allows, in 
principle, to use any functional form for g{x). We note they for the particular case g{x) = x 
e has dimensions of inverse length. To model the laser pulses, we shall use a trapezoidal 



4 



envelope f(t) given by 
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for < t < tl 
for ti<t<t2 
for t2 <t <t3 
elsewhere 



(5) 



where ti = 2'iTnon/u, ^2 = ^1 + 27inp/uj, and ^3 = ^2 + ^irnoff/oj. Hon, and Uoff are the 
number of cycles of turn on, plateau and turn off, respectively. 

The initial state in the ID-TDSE is the ground state (GS) of the system before we turn 
on the laser {t = —00) and it can be found solving an eigenvalue problem once the spatial 
coordinate x has been discretized. The corresponding eigenvalue for the potential ([2]) is 
found to be £gs = —18.2 eV (—0.67 a.u.). For comparison, we note that the ground state 
energy of Ne is -21.6 eV (-0.79 a.u.), and -15.8 eV (-0.58 a.u.) for Ar. 

Equation ([1]) can be solved numerically by using the Crank-Nicolson scheme . In order 
to avoid spurious reflections from the boundaries, at each time step, the total wave function 
is multiplied by a mask function of the form cos^^^, which varies from 1 to starting from 
the 2/3 of the grid Q- 

Once having found the state "^IxA^^oithe system from the ID-TDSE ([T]), we can calculate 
the harmonic spectrum as follows IJ]. The harmonic yield of an atom is proportional to 
the Fourier transform of the acceleration a{t) of its active electron. That it, 
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where a{t) can be obtained by using he following commutator relation 

d2(x) 



a{t) 



i^{t)\['H{t),[H{t),x]]\^{t)), 



(7) 



where 'H(t) is the Hamiltonian defined in the Eq. ([T]). The function D{u) is called the dipole 
spectrum, since D{uj) gives the spectral profile measured in HHG experiments. 



B. The Strong Field Approximation (SFA) for inhomogeneous fields 

Another model to evaluate high-harmonic spectra for atoms in intense laser pulses is the 
Lewenstein model j^. The main ingredient of this approach is the evaluation of the time 
dependent dipole moment d[t). Within the Single Active Electron (SAE) approximation 
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and considering the harmonic radiation is directed mainly in the x-axis d{t) in the length 



gauge for an atom can be written jo] 



4(t) = -i /" dt' [ dkd,onA^ + Mt'),t') 

J to J 

X (i;,^^(k + A(t))exp[-i^o(k,t,t')] +C.C (8) 



5o(k, t, t') = I dt" { ^"""7^ + /, } (9) 



In dHD 

^„l[k + A{t")] 
2 

is the semi classical action, Ip is the ionization potential of the atom, and A{t) = 
— J^^E(t')dt is the vector potential of the linearly polarized laser electric field E(t). The 
ionization and recombination matrices are given by 

ciion,x(k,t) = {^l,\E{t)x\(f)o) (10) 

and 

c?rec,x(k) = (^k| - a;|0o) (11) 

respectively. \E'k is a normalized plane wave of momentum k 

*k(r) = (27r)-=^/V" (12) 

and 00 is the undressed initial-state of the atom. In our studies we use hydrogenic Is states 
of the form 

Mr) = y^e-^^ (13) 

and we chose the effective charge A in order to match the energy of the ground state Sgs 
of the ID model atom (see Sec. II. A), i.e. A = a/2|£^gs|- Using ( !T2|) and (fT3!) the explicit 
expressions for dion,x(M) and drec,x(M) are 

27/2_)^5/2 ^ 

d.onAKt) = i ^ (k^ + A2)3 ^(^) (14) 



and 



27/2_)^5/2 ^ 

drecA^) = , (1,2+V)3 ^^^^ 
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respectively. The spectrum of the emitted hght polarized along the x-axis is obtained by 
modulus squaring the Fourier transform of the dipole acceleration 

a^(fi) = / dtd^{t) exp{int) (16) 
Jo 

where the integration is carried out over the duration of the laser pulse, Tp, by applying 



a fast Fourier transform algorithm. The numerical calculation o 



tidimensional integration over momentum and time. As usual jo], we have performed the 
three-dimensional integration over k using the saddle point or stationary phase method, 
meanwhile all time integrations are performed numerically. 

Eq. ([H]) has a direct interpretation in terms of the three step or simple man's model j^. 
The first step is the strong-field ionization of the atom or molecule as a consequence of 
the nonperturbative interaction with the coherent electromagnetic radiation. The classical 
propagation of the electron in the field defines the second step of the model. Finally, the 
third step in the sequence occurs when the electron is steered back in the linearly polarized 
field to its origin, recombining under the emission of a high-energy photon. One of the main 
features of the HHG process is the coherence of the emitted radiation, which, e.g., opens 



Eq. ([8]) involves a mul- 
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the possibility of generating attosecond pulses 

The Lewenstein model implicitly considers homogeneous laser electric fields, i.e. the elec- 
tric field does not change in the region where the electron dynamics takes place. In order to 
consider non- homogeneous fields, the SFA needs to be modified accordingly. Our goal is to 
find the modifications in the electron momentum and the classical action produced by the 
non- homogeneous field. The first step is to find the electron trajectories starting from classi- 
cal arguments employing the Newton equation for an electron moving in a non-homogeneous 
electric field E{x,t), linearly polarized in the x-axis. It can be then written as 

= -E{t){l + 2ex{t)) (17) 

where Viaser{x, t) is given by ([3]) and we have collected the time dependent part of the electric 
field in E{t), i.e E{t) = Eof{t) sin cut. We use the Picard iteration |l^ extended to second 
order ordinary differential equations to solve Eq. (fT7|) and we restrict ourselves only to the 
first order term. The {n + l)th order solution can be written in terms of the nth one as 



follows 



Xn+l{t) = Xo + Vo{t- to) 



(18) 



to 



fit",X4t"),Xnit"))dt' 



to 



dt' 



where xo = x(to), vq = xo = x(to) and in our case f(t,Xn(t),Xn(t)) = —E(t){l + 2eXn(t)), 
E(t) being the time dependent part of E{x,t), i.e. E(t) = £'o/(t) sinwt. Considering the 
initial conditions, i.e. Xq = and Vq = (the electron starts its movement at the origin with 
zero velocity), we finally obtain 



Xi(t) = a{t)-a{to)-A{to){t-to) 



(19) 



where we have defined a{t) = dt'A(t'). The next step is to calculate the classical action 
and the saddle point electron momentum starting from the electron trajectories. After 
elementary algebra we can write for the modified action S{\i.,t,t') 



S{k,t,t') = ^o(k,t,t') 



+ 2e 



k- / dt"A{t")Xi{t") 
Jt' 

f dt"A\t")X,{t") 
Jt' 



(20) 



where So(k,t,t') is defined in ([9]) and Xi(t) is the electron trajectory of Eq. 
saddle point momentum is found from the stationary condition 



Vk^(k,t,t') = 



H). The 



(21) 



Consequently 



with 



Vk5(k,t,t') = Vk^o(k,t,t') + 2e [ dt"A{t")Xi{t") 

Jt' 



Vk^o(k, t, t') = k(t - t') + a{t) - a{t') 

Finally we obtain for the saddle point or stationary momentum kst(t, t') 

a{t) - a{t') 



(22) 



(23) 



Kt{t,t') 



it - 1') 

2e 
it - f) 



dt"A{t")X,{t") 



(24) 



Replacing ([M]) in ([20]) results 



S{K,t,t') = So{K,t,t') 



2e 



ksf f dt"A{t")X^{t") 
Jt' 

dt" A\t")Xi{t") 



(25) 



The time dependent dipole moment, Eq. (j8]), can be written, after the saddle point method 
for the momentum k is applied jol, as 

3/2 



4(t) 



.i / dt' ( / 



X d,o„,,.(k,t(t,t') + A(t'),t') 

X ci;,,^.(k,i(t,t') + A(t))exp[-i5(k,t,t,t')] +C.C 



(26) 



where r/ is a small parameter. We recover the homogeneous dipole moment putting e = in 

m &• 



III. RESULTS 



Results using the ID-TDSE model developed in Sec. II. A are presented in Figs. 1-7. In 
Figs. 1-6 the laser intensity is / = 2 x 10^^ W-cm~^ and the laser wavelength is A = 800 nm. 
We have used a trapezoidal shaped pulse with two optical cycles turn on, i.e. Uon = 2, and 
turn off, i.e. Uoff = 2, and a plateau with six optical cycles, i.e. rip = 6 (10 optical cycles in 
total, i.e. approximately 27 fs). The model atom has Scs = —0.67 a.u. and it is driven by 
a laser electric field of the form E{x,t) = £'o/(t)(l +ex(t)), where Eq is the peak amplitude 
and f{t) is the pulse envelope (see panel (e) of Fig. 2 for a plot of the laser electric field). 

One additional parameter we take into account in the ID-TDSE simulations is the spatial 
region where the electron moves. This parameter will naturally appear in real situations 
using bow-tie- shaped nanostructures, as those employed in the experiments of Kim et.al 
and when linearly polarized pulses are utilized, that restricts the electron dynamics mostly 
to one dimension. The bow-tie shaped nanostructures are characterized by a spatial gap that 
can be, in principle, changed between certain ranges when the nanostructure is engineered. 
We use in our ID-TDSE model spatial grids in terms of the quiver radius, = Eq/u"^, u 
being the driven laser frequency, {u = 0.057 a.u. for a laser wavelength A = 800 nm) that in 
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the case under study is ao ~ 23.2 a.u. (1.23 nm). Three different spatial grid sizes will be 
employed, namely xum = ±7. Sao, = ±4.5ao and xum = il.Sao, corresponding to gaps 
of 18.7 nm, 11 nm and 3.7 nm, respectively . 

In Fig. 1 we plot the HHG spectra for different values of e and for a grid size of xum = 
±7.5ao- Panel (a) is the homogeneous case and we have chosen values of 0.01 (panel b), 
0.02 (panel c) and 0.05 (panel d) for the inhomogeneity parameter e that corresponds to 
inhomogeneity regions of 100 a.u. (5.3 nm), 50 a.u. (2.7 nm) and 20 a.u. (1 nm), respectively. 
Two distinct characteristics can be observed in panels (b)-(d), and we can summarize them 
as follows (i) Odd and even harmonics are now present. This new feature appears because we 
have broken the symmetry of the total potential, Vatom + Viaser, introducing a new asymmetric 
term, i.e. ex'^; (ii) The absence of a clear HHG cutoff. This effect is related with the electron 
trajectories that contribute to the harmonic spectra and will be object of study using the 
Gabor analysis and the semiclassical simulations (see below for details). 

Our next step is to use time-analysis tools in order to characterize the HHG spectra 
calculated using the ID-TDSE model. To this end we employ the Gabor transformation, 
developed in the 1940s by D. Gabor [l8|, that has shown to be a very powerful tool in order 
to estimate the emission times of HHG in atoms and molecules and to discriminate the 
different electron trajectories [l^]. Starting from the dipole acceleration a{t) of Eq. ([7]) the 
Gabor transform is defined as 

aGin,t) = [ dfaif f""^ ^"^^ - tOV2^'] e^p(,^^/) (27) 
J crv27r 

where the integration is usually taken over the pulse duration. In our studies we use a = 
l/3a;, with u being the central laser frequency. With this value of a we achieve an adequately 
balance between the time and frequency resolutions (see Ref. 17|] for details). Results of 
the Gabor analysis of the HHG spectra of Fig. 1 are presented in Fig. 2. From panels 
(a)-(d) can be observed the clear differences between the homogeneous (panel (a)) and non- 
homogeneous cases (panel (b)-(d)). In the zoomed regions we show a time interval during 
the laser pulse for which a complete electron trajectory, from birth time to recollision time, 
falls within the pulse plateau. From these plots it is possible to clearly observe the short 
(for emission times smaller than ~320 a.u.) and long (the emission times are larger than 
~320 a.u. for this case) trajectories (see Ref. [l7| for more details). The highest harmonic 
order (around the 40th) is generated at an emission time of ~320 a.u. for the homogeneous 
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case. On the other hand, only an extended short trajectory (with emission times smaller 
than ~320 a.u.) is present for all the nonhomogeneous cases and no a clear harmonic cutoff 
is visible. 

Alternatively, in Figs. 3 and 4 we plot the HHG spectra for different values of e and 
the Gabor analysis, respectively. Here a grid size of xum. = ±4.5ao is used. Panel (a) is 
the homogeneous case and we have chosen values of 0.01 (panel b), 0.02 (panel c) and 0.05 
(panel d) for the inhomogeneity parameter e that corresponds to inhomogeneity regions of 
100 a.u. (5.3 nm), 50 a.u. (2.7 nm) and 20 a.u. (1 nm), respectively. We note that the 
difference between these last two figures (Figs. 3 and 4) and Figs. 1 and 2 is hardly visible 
for all the cases, showing that the decrease of the grid size has no effect in the HHG spectra, 
both for the homogeneous and nonhomogeneous cases. 

Finally, in Figs. 5 and 6 we plot the HHG spectra for different values of e and the 
Gabor analysis, respectively, and now with a grid size of xum = ±1.5ao. Panel (a) is the 
homogeneous case and we have chosen values of 0.01 (panel b), 0.02 (panel c) and 0.05 
(panel d) for the inhomogeneity parameter e that corresponds to inhomogeneity regions of 
100 a.u. (5.3 nm), 50 a.u. (2.7 nm) and 20 a.u. (1 nm), respectively. Here it is possible to 
observe the substantial difference between panels (b)-(d) of these last two figures and those 
from Figs. 1-4, showing that the electron confinement, joint with the inhomogeneities of the 
laser electric field, are the responsible of the appearance of a clear extended harmonic cutoff. 
We also note the HHG spectra for the homogeneous case (panel (a) of Figs. 1, 3 and 5) are 
identical for all the cases, confirming that not only the electron confinement is the reason of 
the new features present in the HHG spectra. In some sense, when we restrict the electron 
movement, only the short trajectories contribute to the formation of the harmonic spectrum 
and this feature is clearly visible in the zoomed regions of the Gabor analysis in Fig. 6 where 
only short trajectories are distinguishable, i.e. those with emission times smaller than ~320 
a.u. 

In order to explore how our ID-TDSE approach behaves we have calculated HHG spectra 
for longer wavelengths. We present our results in Fig. 7 for A = 3.2 /im. This value for A 



represents an example that could be perfectly reachable experimentally |19l . |20| and allow 
us to reach values for the electron excursion region (i.e. the quiver radius Oq) closer to 
those considered in real experiments The laser intensity is / = 1 x 10^^ W-cm~^ and we 
have employed a gaussian shaped pulse with 6 cycles FWHM, f{t) = exp , where 
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T is the FWHM (full-width at half-maximum) duration of the driving laser intensity I{t) 
(oc l-E(t)p). We use a spatial grid of xum = ±1.5ao, where ao = 83 a.u. (4.4 nm), i.e. 
the electron dynamics is confined in a region of around 13.2 nm. The panels correspond to 
the homogeneous case (a), e = 0.01 (b), e = 0.02 (c) and e = 0.05 (d), respectively. We 
have included zoomed panels for the (a) and (d) cases. These two plots allow us to observe 
clearly the presence of odd and even harmonics for this particular nonhomogeneous case and 
only the appearance of odd harmonics for the homogeneous one. From Fig. 7 we conclude 
that for longer wavelengths the cutoff extension due to the nonhomogeneities of the field, 
combined with the electron confinement, is far more pronounced, e.g. for the case of A = 3.2 
fim and e = 0.05 the cutoff is almost 3 times larger than the homogeneous case. This 
extension could open the avenue for the production of attosecond pulses in the keV regime 
(for a theoretical study at A = 800 nm see llj). Furthermore the region of confinement 
using these laser parameters is closer to values that could be found in real nanostructures. 

In the following we use the modified SFA model presented in the Sec. II. B to generate 
HHG spectra produced by nonhomogeneous fields. We employ a laser pulse with an electric 
field of the form E{x,t) = Eof{t){l + ex(t)), where Eo is the peak amplitude, /(t) is the 
pulse envelope and we use only the first order for x{t), i.e. x{t) ~ Xi{t) with Xi{t) of (fT9l) . 
In Fig. 8 we show the predictions of this model and Fig. 9 represents the Gabor analysis 
of the HHG spectra of the former. In order to enhance the new HHG features and see 
more clearly the differences between the studied cases now the laser field has an intensity 
/ = 6 X 10^^ W-cm~^ and a wavelength of A = 800 nm. We have used a trapezoidal shaped 
pulse with three optical cycles turn-on [uon = 3) and turn off {uoff = 4) and a plateau 
of four constant-amplitude optical cycles {up = 4) (see panel (e) of Fig. 9) and the model 
atom has Sgs = —0.67 a.u. With these parameters the predicted cutoff is Uc = 85 as it is 
shown by the arrow plotted in panel (a). Panels (b), (c) and (d) correspond to values of 
2e = 0.01, 2e = 0.02 and 2e = 0.05, respectively. As in the case of the ID-TDSE model, the 
SFA approach predicts the appearance of odd and even harmonics and a cutoff extension. 
In order to show more clearly these features we have zoomed out a small region of the 
HHG near the cutoff for panels (a) and (b). In this graph it is more transparent to see 
how the inhomogeneous fields modify the harmonic spectra. An additional feature appears 
and it could be extracted from panels (a) (homogeneous case) and (d) {2e = 0.05): there 
exists a clear enhancement in the harmonic yield in the plateau region (see e.g. a region 
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around the 25th-40th harmonics). The Gabor analysis of the HHG spectra calculated using 
the SEA, Fig. 9, shows similarities and differences between this model and the ID-TDSE 
approach (the SFA approach does not account for the restrictions on the electron motion). 
For example, in panel (c) only short trajectories appear to contribute to the HHG spectra 
(similar to the xum = ±1.5ao case), meanwhile panel (d) has analogous features to panels 
(d) of Figs. 2 and 4 (i.e. the xum = ±7. Sao and xum = ±4.5ao cases respectively). This 
behavior could be related to the limited validity of the perturbation-like character of the SFA 
approach developed in Sec. 11. B. An extension of the SFA approach including the quantum 
orbits analysis will be published elsewhere. 

We now concentrate our efforts in to explain, using the semiclassical three-step model, 
the extension of the harmonic cutoff. As was already pointed out, this new feature appears 
as a consequence of the combination two factors, namely the nonhomogeneous character of 
the laser electric field and the electron confinement. It is well known that the position of 
the HHG cutoff holds 

ricU = 3.17Up + Ip (28) 

where ric is the harmonic order at the cutoff, u the laser frequency. Up the ponderomotive 
energy [Up = I/Au"^, I being the laser intensity in a.u.) and Ip the ionization potential of 
the atom or molecule 9|]. This relationship can be obtained solving the classical Newton 
equation for an electron moving in a linearly polarized electric oscillating field under the 
following conditions (the resulting model is also known as the simple man's model): (i) the 
electron starts with zero velocity at position zero at time t = to (to is known as ionization 
time), i.e. x{to) = and x{to) = for our ID model; (ii) when the electric field reverses its 
direction, the electron returns to its initial position (i.e recollides with its parent ion) at a 
later time t = ti (ti is also known as recollision time), i.e. x{ti) = 0. The electron kinetic 
energy at the return time ti can be obtained from Ek{ti) = x{tif' /2 and finding the value 
of ti (as a function of to) that maximizes this energy, Eq. (l28l) is fulfilled. 

We have solved numerically the Newton equation for an electron moving in a linearly 
polarized (in the x-axis) electric field with the same parameters used in the ID-TDSE 
calculations. Specifically we find the numerical solution of x{t) = —'VxViaseriXjt) with 
Viaser{x,t) givcu by ([3]) and E{x,t) of Eqs. (jlj) and ([5]) without any approximation, i.e. 
we solve the Eq. (fT7|) . For fixed values of ionization time to it is possible to compute the 
classical trajectories and to numerically calculate the times ti where the electron recollides, 

13 



i.e. x(ti) = 0. Also, once the ionization time to is fixed, the electron trajectory is completely 
determined. In Fig. 10 panels (a)-(d) we plot the dependence of the harmonic order on 
the ionization time (to) and recoUision time (ti), calculated from n = {Ek{ti) + Ip)/uj, with 
i = and i = 1, and for the cases presented in the ID-TDSE simulations, i.e. homogeneous 
(panel a), e = 0.01 (panel b), e = 0.02 (panel c) and e = 0.05 (panel d), respectively. At 
this point we have not restricted the electron trajectories, and consequently we allow the 
electron to move in all the space (see below for details). The temporal axis, i.e. the x-axis, 
is plotted in terms of optical cycles and we have chosen a temporal window from 3.5 to 5 
optical cycles (i.e. from 380 a.u. to 550 a.u.). In all the panels both the short and long 



trajectories (see see e.g. 



Ill) are labeled. From panel (a) it is possible to observe that the 



maximum kinetic energy of the returning electron is in perfect agreement with Eq. 
(no harmonic order beyond ric ~ 36 is reached). On the other hand panels (b)-(d) show 
how the nonhomogeneous field modifies the electron trajectories and that no clear HHG 
cutoff is observed. This is in consistent with the predictions of the ID-TDSE for the largest 
spatial grids (see Figs. (2) and (4)). Although from panels (b)-(d) we observe a no clear 
HHG cutoff, the kinetic energy gained by the electron is finite in agreement with the energy 
conservation. Furthermore, panels (b)-(d) show similar features as those observed in Fig. (2), 
panels (b)-(d) and Fig. (4), panels (b)-(d), i.e. only extended short trajectories contribute 
to the harmonic radiation. This characteristic is related with the modifications the electron 
trajectories suffer due to the nonhomogeneities of the field (see below for details). 

In order to classically simulate the ID-TDSE results, but for the smallest grid size, 
i.e. xiim = ±1.5ao, we restrict the classical electron trajectories to the domain [— aojC^o]- 
The ±ao values represent the starting point of the mask function and consequently a fair 
comparison is possible. The results are presented in Fig. 11, panels (a)-(d). From these 
plots we can argue that only short trajectories contributes to the harmonic radiation. This 
is related with the electron motion restriction, i.e. the confinement, we have incorporated 
in the classical simulations. Furthermore, a clear HHG cutoff is now observed for all the 
nonhomogenous cases and its value is in clear agreement with the ID-TDSE predictions (see 
Figs. (5) and (6)). 

Finally, in Fig. 12 the recollision time ti of the electron is presented as a function of 
the ionization time to and for several values of e. Panel (a) represents the non-confined 
case and in panel (b) we have restricted the electron motion into the region [— ao,ao]- The 
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long trajectories are those with recolhsion times ti > 4.25 optical cycles and only for the 
homogeneous case (red squares (□)) these trajectories are clearly visible. On the other hand, 
short trajectories are characterized by ti < 4.25 optical cycles and these are present for both 
the homogeneous and nonhomogeneous cases. Our results are consistent with those shown 
in Uj, but note our inhomogeneity parameter is more than one order of magnitude higher. 
From panel (a) we observe how the long trajectories are modified by the nonhomogeneity, 
namely the homogeneous long trajectories (red squares (□)) with ionization times to around 
the 3.25 and 3.75 optical cycles merge into unique trajectories. Additionally the branch with 
to ~ 3.75 has now ionization times more than half an optical cycle smaller when e increases; 
hence, the times spent by the electron in the continuum increase. This fact explains the 
vanishing long trajectories seen in the panels (b)-(d) of Fig. 10. Additionally, the electric 
field strength at the ionization time for short trajectories is higher than for long trajectories 



and considering the ionization rate as a highly nonlinear function of this electric field [2lU22|. 
long trajectories are much less efficient than the short ones. On the other hand, short 
trajectories are almost independent of e and only for higher values noticeable differences are 
visible. When the electron motion is confined, panel (b), only short trajectories are present 
for all the cases and this confirm the fact that long trajectories are absent and only the short 
trajectories contribute to the formation of the harmonic radiation. 



IV. CONCLUSIONS AND OUTLOOK 



We studied high-order harmonic generation in a model atom produced by non- 
homogeneous fields. These fields are not merely a theoretical speculation but are present 
in a vicinity of a metal nanostructure when it is irradiated by a short laser pulse. We have 
extended two of the most widely used models, namely the numerical solution of the time 
dependent Schrodinger equation (TDSE) in reduced dimensions and the semiclassical ap- 
proach known as Strong Field Approximation (SFA). Both models are able to predict the 
new features that appear due to presence of inhomogeneities, namely the appearance of odd 
and even harmonics and an extension in the cutoff position. The latter feature is a conse- 
quence of the combination of the nonhomogeneous electric field and the confinement of the 
electron dynamics as can be extracted from the Gabor analysis and the semiclassical simula- 
tions made for all the studied cases. We have proved the robustness of the ID-TDSE model 
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using different laser parameters and solving the classical equation of motion for an electron 
moving in a linearly polarized nonhomogeneous electric field we are able to understand the 
reasons of the cutoff extension. 

The models presented allow to use, in principle, any functional form for the nonhomoge- 
neous fields. Moreover they are able to incorporate real parameters in the simulations, such 
as the region dimensions where the electron dynamics takes place. We plan to continue our 
investigations in this direction, joint with studies of other laser-matter processes driven now 
by these nonhomogeneous fields. 
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Figures captions 



Fig. 1. High-order harmonic generation (HHG) spectra for a model atom with Eqs = 
—0.67 a.u. generated using the ID-TDSE model and with a spatial grid of xum = ±7.5ao 
(see text for details). The laser parameters are / = 2 x 10^^ W-cm~^ and A = 800 nm. We 
have used a trapezoidal shaped pulse with two optical cycles turn on, i.e. Hon = 2, and turn 
off, i.e. Uoff = 2, and a plateau with six optical cycles, i.e. Up = 6 (10 optical cycles in 
total, i.e. approximately 27 fs). The arrow indicates the cutoff predicted by the semiclassical 
model jsl. Panel (a) homogeneous case, (b) e = 0.01 (100 a.u), (c) e = 0.02 (50 a.u) and (d) 
e = 0.05 (20 a.u). 

Fig. 2. (Color online) Panels (a)-(d): Gabor analysis for the HHG spectra of Figure 
1. The zoomed regions in all panels show a time interval during the laser pulse for which 
the complete electron trajectory, from birth time to recollision time, falls within the pulse 
plateau (see Ref. [17| for details); panel (e) shape of the laser electric field. In panels (a)-(d) 
the color scale is logarithmic. 

Fig. 3. Idem Fig. 1 but with a spatial grid of xum = ±4.5ao. 

Fig. 4. (Color online) Gabor analysis for the HHG spectra of Figure 3. The zoomed 
regions in all panels show a time interval during the laser pulse for which the complete 
electron trajectory, from birth tinre o recolli=io,r time, talis wrthi,r the pulse plateau (see 
text and Ref. [17| for details). 

Fig. 5. Idem Fig. 1 but with a spatial grid of xum = ±1.5ao- 

Fig. 6. (Color online) Gabor analysis for the HHG spectra of Figure 5. The zoomed 
regions in all panels show a time interval during the laser pulse for which the complete 
electron trajectory, from birth time o recollision time, falls within the pulse plateau (see 
text and Ref. jl3] for details). 

Fig. 7. HHG spectra for a model atom with Sgs = —0.67 a.u. generated using the 
ID-TDSE model and with a spatial grid of xum = ±1.5ao (see text for details). The 
laser parameters are / = 1 x 10^^ W-cm~^ and A = 3.2 fim (3200 nm). We have used a 
gaussian shaped pulse with 6 cycles FWHM. The arrow indicates the cutoff predicted by 
the semiclassical model j9|. The insets in panels (a) and (d) show particular zoomed regions 
of the harmonic spectra near the cutoff region (see the text for details). 
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Fig. 8. (Color online) HHG spectra for a model atom with £gs = —0.67 generated using 
the extended SFA approach. The laser parameters are J = 6 x 10'^'^ W-cm~^ and A = 800 
nm. We have used a trapezoidal shaped pulse with 3 optical cycles 'turn-on' and 'turn-off' 
and a plateau of 4 const ant- amplitude optical cycles (see panel (d) of Fig. 9). The arrow 
indicates the cutoff predicted by the semiclassical model |9|. Panel (a) homogeneous case, 
(b) 2e = 0.01 (100 a.u), (c) 2e = 0.02 (50 a.u) and (d) 2e = 0.05 (20 a.u). The zoomed 
regions correspond to harmonic order values between 100 u and 120 u in panel (a) and 
330 u and 350 u in panel (d), respectively (see text for further details). 

Fig. 9. (Color online) Panels (a)-(d): Gabor analysis for the HHG spectra of Figure 8. 
The zoomed regions in all panels show a time interval during the laser pulse for which 
the complete electron trajectory, from birth time to recoUision time, falls within the pulse 
plateau (see text and Ref. [l^ for details); panel (e) shape of the laser electric field. In 
panels (a)-(d) the color scale is logarithmic. 

Fig. 10. (Color online) Total energy of the free electron (in terms of the harmonic order) 
in the laser field when it recoUides with its parent ion obtained from Newton's second law 
and plotted as a function of the ionization time (green open circles) or the recollision time 
(red filled circles). Panel (a) homogeneous case, (b) e = 0.01 (100 a.u), (c) e = 0.02 (50 a.u) 
and (d) e = 0.05 (20 a.u). In all the cases the motion of the electron is not restricted. 

Fig. 11. (Color online) Idem Fig. 10, but with the motion of the electron confined into 
a region [— ao, +ao]- 

Fig. 12. (Color online) Dependence of the semiclassical trajectories on the ionization 
and recollision times for different values of e and for the non confined case, panel (a) and 
the confined case, panel (b). Red squares (□), homogeneous case e = 0; green circles (Q) 
e = 0.01; blue triangles (A) e = 0.02 and magenta diamonds (o) e = 0.05. 
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